Fourier law in the alternate mass hard-core potential chain 
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We study energy transport in a one-dimensional model of elastically colliding particles with al- 
ternate masses m and M. In order to prevent total momentum conservation we confine particles 
with mass M inside a cell of finite size. We provide convincing numerical evidence for the validity of 
Fourier law of heat conduction in spite of the lack of exponential dynamical instability. Comparison 
with previous results on similar models shows the relevance of the role played by total momentum 
conservation. 
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After several decades of intensive investigations [1-9], 
the precise conditions that a dynamical system of inter- 
acting particles in ID must satisfy in order to obey the 
Fourier law of heat conduction are still not known. 

For non-interacting particles in external potential it 
has been shown [3] that exponential local instability leads 
to Fourier law. Actually, even linear mixing without ex- 
ponential instability, such as found in generic polygonal 
billiards [10], has been shown to be sufficient for a dif- 
fusive heat transport [4] . In addition, several interacting 
non-integrable many-particle systems which clearly obey 
the Fourier law have been proposed and investigated [1] . 
However, it should be noted that in the above models, 
the total momentum is not conserved. In several recent 
papers [5,6,8] it has been suggested that total momentum 
conservation does not allow Fourier law. Moreover, us- 
ing renormalization group [8] , it is argued that a generic 
momentum conserving particle chain should, in a macro- 
scopic limit, be equivalent to ID hydrodynamics with 
thermal noise where the coefficient of thermal conductiv- 
ity should diverge with the system size L as k(L) oc L^^^. 
However, most existing numerical data do not support 
this universal constant. Instead, it has been proposed 
that k{L) oc L2-2/q j-gj^ where a is the exponent of the 
diffusion ((Aj;2) 2Dt"), <a<2. 

We would like to remark that there exist a model 
[11] (a particle chain with interparticle potential V{x) = 
1 — cos(a::)) in which, in spite of momentum conservation, 
the heat conduction seems to obey the Fourier law. The 
reason for such behaviour is not clear and the precise role 
of the total momentum conservation needs to be clarified. 

In previous papers two models have been considered, 
both are mixing and without exponential instability: (i) 
the triangular billiard channel [4] , which exhibits Fourier 
law and (ii) the alternate mass hard point gas model [7] 
in which the coefficient of thermal conductivity diverges 



with the system size. The difference between the two 
models is that in case (ii) the total momentum is con- 
served while in case (i) it is not. 

On the other hand, model (i) has been criticized since, 
in spite of the fact that one can perfectly well define an 
internal local temperature, there is no mechanism to pro- 
vide local thermal equilibrium due to lack of interparticle 
interaction, and this may look somehow unsatisfactory. 
In this paper we consider a model which is identical to 
the alternate mass hard-point gas [7], namely it consists 
of a one-dimensional chain of elastically colliding parti- 
cles with alternate masses m and M . Here, however, in 
order to prevent total momentum conservation we con- 
fine the motion of particles of mass M (bars) inside unit 
cells of size I. Schematically the model is shown in Fig. 
1 in which particles with mass m move horizontally and 
collide with bars of mass M which, besides suffering col- 
lisions with the particles, are elastically reflected back 
at the edges of their cells. In between collisions, parti- 
cles and bars move freely. Our numerical results clearly 
indicate that our model, contrary to the translationally 
invariant model [7], obeys the Fourier law. The only 
difference between the two models is total momentum 
conservation. 
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FIG. 1. The geometry of the model. The elementary 
cell (indicated by two dotted lines) has unit length 1 = 1. 
The bars have mass M = 1, and the particles have mass 
m = (\/5 — l)/2. The two heat baths at temperatures Tl and 
Tr are indicated. 
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The total length of the system is L = Nl where N is 
the number of fundamental cells. In all the calculations 
presented in this paper, we fix Z = 1 so that L — N. We 
also take M = 1, m = {V^— l)/2 and we verify that the 
numerical value of the mass ratio is not relevant. 

A direct way to test whether the system obeys the 
Fourier law is to put two heat baths with small temper- 
ature difference into contact with the two ends of the 
system, and check the dependence of the thermal cur- 
rent on the system size. Here statistical thermal baths 
are used; that is when the first (last) bar collide with the 
left (right) side of the first (last) cell, it is injected back 
with a new speed generated from the distribution 

M\v\ ( Mv^ \ 
Pl,r{^) = exp -— — . (1) 

This assures, for small temperature gradients, that the 
edge particles have canonical (Maxwellian) velocity dis- 
tribution. In our simulations we fixed Tl — 1.1 and 
Tr = 0.9. 
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FIG. 2. Internal local temperature as a function of the 
rescaled position x/N. The total number A'' of cells is: A'^ = 64 
(dotted line), N = 128 (solid line), and = 256 (dashed 
line). Notice the good scaling behavior of the temperature 
field. The inset shows kurtosis of the local velocity distribu- 
tion. 



the average time measured by the total collisions number 
is larger than 5 x 10^° for N = 512 which is the largest 
system size we have considered. 

We should emphasize that in previous models [3,4,7] 
the local thermal equilibrium cannot be established, 
whereas in the model considered here, the local ther- 
mal equilibrium is well established independently of the 
thermal baths used. We have checked that the velocity 
distribution function for each bar and particle is a Gaus- 
sian function whose width gives the local temperature, 

whereas the kurtosis, defined by K{x)^^'') = ^^i^2 — 1, 

for bars and K{x^^'^^) = ^j^y — 1, for particles, are close 

to zero. The kurtosis versus the bar/particle site is plot- 
ted in the inset of Fig. 2. 

In Fig. 3 we show the stationary time-averaged heat 
fiux (j) as a function of the system size N. The best fit of 
numerical data gives (j) = 0.24Af-'' with 7 = 0.99±0.01. 
The coefficient of thermal conductivity appears there- 
fore to be independent on TV, which means that the 
Fourier law is obeyed. Its numerical value reads as 
K = = 1.20 ±0.05. 




FIG. 3. Scaling behavior of the stationary time-averaged 
heat fiux (j) as a function of the system size A*'. The least 
square fit gives a slope —0.99 ± 0.01. 



For any given initial condition, after a long enough 
transient time, the system reaches a stationary state. 
Then one may compute the local temperature, defined 
as 

T{x) = (Mv'^) for X = x^" = {i - 0.5), i = l,---N 
T{x) = (tou^) for X = a;P"* = j, j^l^-.-N. (2) 

Here x\^'' and x^^"^* can be regarded roughly as the time- 
averaged positions of the ith bar and the jth particle. 

In Fig. 2 we plot the temperature profile versus the 
scaled length x/N for different values of N . Notice the 
good hnear scaling behaviour. In both Fig. 2 and Fig. 3, 



To investigate how the energy diffuses along the sys- 
tem, we set Tl = Tr ~ 1. Then, after the equilibrium 
state is reached, the middle particle (i.e. the 100th one 
for N = 200 in Fig. 4a) is given a speed u = 1/^/ra, so 
that its energy is four times bigger than the equilibrium 
average. The evolution of the energy profile along the 
chain is then recorded afterwards. To suppress statisti- 
cal fluctuations, 10^ realizations are taken into account 
for the average. The width of the energy profile can be 
measured by its second moment 

2 _ l{E{x,t)- Ea)(x-xafdx 
^^'^ J{Eix,t)-Eo)dx ' ^' 
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where E{x,t) = *)+^, Mfsix-x^^"^). 

In our calculations for Fig. 4, Eq — 0.5 and xq = 100. 
The energy profile spreads as a'^{t) = 2Dt with D = 
1.20 ± 0.01 (Fig. 4b), which agrees with the thermal con- 
ductivity K = 1.20 ± 0.05 very well. 
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FIG. 5. The Green-Kubo integral versus t for N — 24 and 
T = 1. To get (J(t) J(0)) used in (6), a single orbit is followed 
up to time tt — 4.8 x 10^; the error, S{J{t)J{0)), is estimated 
in the same way as in Fig. 6. Dotted lines (almost indistin- 
guishable) indicate the error boundaries of n{t), obtained by 
replacing the integrand in (6) by {J(r)J(0)) ± S{J{t)J{0)) 
respectively. The errors are very small and negligible. 
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FIG. 4. Diffusive property of our model. Initially the sys- 
tem is at equilibrium with temperature T = 1; then an im- 
pulse is given to the particle in the middle of the chain, (a) 
Energy distribution along the channel at different times, (b) 
The second moment of the energy distribution versus time. 

If the system obeys the Fourier law, its thermal con- 
ductivity can also be obtained via the Green-Kubo for- 
mula 

K^ \im hm f\j{T)JmdT, (4) 

where the heat current can be written as [12] 

^ 1 

1=1 

In applying Green-Kubo theory, a periodic boundary 
condition is imposed to the system. 



In Fig. 5 we plot the quantity 

<t) ^ j\jiT)J{0))dT, (6) 

as a function of time t ioi N = 24. Since (J(t) J(0)) 
decays in time very fast (see Fig. 6), one has that K{t) 
tends to k very fast as well. The numerical result gives 
K = 1.22, in excellent agreement to the heat conductivity 
obtained via simulations with thermal baths (/c = 1.20). 




FIG. 6. Absolute heat current correlation function for 
N = 2 and T = 1. The time average is taken over a sin- 
gle orbit up to time t; = 5 x 10^. The dotted line shows the 
standard deviation of {J(r) J(0)) estimated over {ti — t)/5to 
values of J(r -I- ro) J(ro) obtained by changing tq with a step 
Stq = 0.05 along the orbit. The dashed line indicates an 
exponential decay exp(— 0.38r). 
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In Fig. 6 we plot the absolute value of current-current 
time correlation function C(t) = (J(r) J(0)) versus r for 
N — 2. It is interesting to remark that numerical re- 
sults seem to indicate a clear exponential decay of corre- 
lation C(t) which seems in contradiction with the linear, 
marginally unstable, dynamics. Indeed, as shown in [10], 
the existence of periodic orbits in a marginally unstable 
system necessarily implies an asymptotic power-law de- 
cay of correlation. The asymptotic power-law tail may 
be however very difficult or impossible to observe numer- 
ically. In fact, what we see is a transient exponential 
decay over several orders of magnitude which is very ro- 
bust against changing system parameters, such as m/M 
or N. 

Finally, we study the behaviour of the thermal con- 
ductivity K versus the mass ratio m/M (Fig. 7). It is 
interesting to note that even at m = M = 1 one finds a 
finite conductivity k, which is close to the maximum of 
the curve in Fig. 7, in spite of the fact that in this case 
(m = M) the dynamics is pseudo-integrable since, for 
the isolated system, the set of magnitudes of the veloci- 
ties of initial particles is conserved. However, the system 
m = M is not strictly integrable since the topology of 
invariant surfaces is more complex than the one of the 
tori. Notice also that only local thermal equilibrium is 
absent for m = M. 




FIG. 7. The thermal conductivity k (bullet) versus the 
mass m for N — 24 and T — 1. AI — 1. The solid curve 
is drawn to guide the eyes. 

In the present paper we have demonstrated diffusive 
energy transport and Fourier law for a marginally sta- 
ble (non-chaotic) interacting many-particle system. We 
have thus clearly demonstrated that exponential instabil- 
ity (Lyapunov chaos) is not necessary for the establish- 
ment of the Fourier law. Furthermore, we have shown 
that breaking the total momentum conservation is cru- 
cial for the validity of Fourier law while, somehow sur- 



prisingly, a less important role seems to be played by the 
degree of dynamical chaos. 
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